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We consider ramifications of the use of high speed light modulators to questions of correlation 
and measurement of time-energy entangled photons. Using phase modulators, we find that temporal 
modulation of one photon of an entangled pair, as measured by correlation in the frequency domain, 
may be negated or enhanced by modulation of the second photon. Using amplitude modulators we 
describe a Fourier technique for measurement of biphoton wave functions with slow detectors. 
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It is the intent of this Letter to explore the ramifica- 
tions of the use of high speed light modulators to ques- 
tions of correlation and measurement of time-energy en- 
tangled photons We start by considering the sys- 
tem of Fig. [T] where a monochromatic pump generates 
non-degenerate time-energy entangled photon pairs. The 
signal and idler photons pass through sinusoidal phase 
modulators. These modulators are driven at the same 
modulation frequency and are connected by a cable such 
that their relative phase may be varied. After passing 
through the modulators the signal and idler photons are 
dispersed, for example by a prism, and the relative posi- 
tions of the signal and idler photons are correlated. When 
the modulation frequency is small as compared to the 
spectral bandwidth of the signal or idler, we find a con- 
sequence of time energy entanglement that we term as 
nonlocal modulation. Specifically, the modulators act cu- 
mulatively to determine the apparent modulation depth. 
For ideal phase modulators the depth of modulation is 
the sum of the modulation depths of the signal and idler 
modulators. When the modulators are run with the same 
phase, the modulation depths add, when they are run in 
phase opposition, the modulation depths subtract; two 
distant modulators with the same modulation depth and 



phase modulators 



pump 




opposite phase will have the same (delta function) fre- 
quency correlation as when both modulators are absent. 

In the latter portion of this Letter, (Fig. 5), we con- 
sider synchronously driven amplitude modulators in both 
channels. We show how they may be used to measure the 
waveform of biphotons that are too short to be measured 
by present day photo-detectors, but are not so short as 
to be out of the range of high speed (60 GHz) light mod- 
ulators. 

Nonlocal modulation as described here may be thought 
of as the time-frequency counterpart of nonlocal disper- 
sion compensation as suggested by Franson Q, and de- 
veloped experimentally by Shih [3| and Silberberg 0, 0] 
and colleagues. Both phenomena depend on the quan- 
tum mechanical addition of probability amplitudes, and 
both do not have a classical analog. There is consider- 
able theoretical and experimental work that is pertinent 
to the work described here. Examples include quantum 
interference effects (f|, nonlocal pulse shaping 0], de- 
scriptions of time-energy entanglement experiments 
showing higher-dimensional entanglement using trains of 
pumping pulses [9|, and experimental control of the joint 



spectrum of down-converted photons [id. 111 | . The wave- 
form measurement technique described below is in part 
motivated by recent demonstrations of EIT based tech- 
niques for generating paired photons that are many cycles 
long at the modulation frequency 12|, |l3|, [14 1 . 

We develop the theory in the Heisenberg picture where 
the output of the parametric down-converter is described 
by frequency domain operators a s (ui) and a](u>j), with 
uji = L0 p — uj, where uj p is the frequency of the monochro- 
matic pump. These operators are expressed in terms of 
the vacuum fields at the input of the nonlinear generat- 
ing crystal a 8 (u, 0) and a](Wj, 0) by frequency dependent 
functions A(ui), B(uj),C(ui) and D(ui), all expressed in 
terms of the signal frequency w. Thus 



FIG. 1: (color online) . Quantum modulation. Signal and idler 
photons are phase modulated at the same frequency and with 
controllable phase. The signal and idler beams are diffracted 
to linear arrays and the positions of the detected photons on 
the photodetectors are correlated. 



a s (oj) =A{uj) a s (oj,0) + B(u) aj(^,0) 
a\(wi) = C{uj) a s (u,Q)+D(u) aj(wj,0). (1) 

The frequency domain operators are the Fourier 
transform of time domain operators, i.e., a(u>) = 
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FIG. 2: (color online). Spectrum of the signal frequency after 
passing through the modulator. In part (a) the modulation 
frequency uu m = 0.1 and is small as compared to the linewidth 
of one unit of the biphoton. In part (b) uo m — 10 and is large 
as compared this linewidth. Counts are per bandwidth in a 
gatewidth T. 

l/(27r) a(t) exp(iuit)dt, that in turn have the com- 
mutator [dj(ti, 0), at(t2, 0)] = S(ti — ti)5j^k- The op- 
erator a(t) is normalized so that the paired count 
rate for a single transverse mode is R = (a<(t)a{t)). 
The coefficients are related by the unitary conditions 
A(w)C*(w) = B(u)D*(lu),\A{uj)\ 2 - \B(u)\ 2 = 1, and 



la . lil have shown 



\D(u>)\ 2 - \C(uj)\ 2 = 1. Gatti et al. 
that for low parametric gain where there is, most often, 
only a single photon at both the signal and idler frequen- 
cies, and with |^4(o;)| = 1, the biphoton is described by 
the function cj)(uj) = A{w)C*{u). 

The phase modulators in the signal and idler channels 
are assumed to be periodic with the same modulation fre- 
quency uj m and with Fourier series ^ n q n exp(— inuj m t), 
and J2 m r ™ exp(— imco m t), respectively. In the frequency 
domain these modulators are described by 

oo 

m s (uj) = ^ q n 6 {w - riL} m ) 

n— — oo 
oo 

m t (uj) = r m S (w - muj m ) . (2) 



Time varying modulators multiply the incoming tem- 
poral waveform, and their transforms convolve with the 
Fourier transform of this waveform. The operators at the 
output of the signal and idler modulators are obtained by 
combining Eqs. (frj and Eqs. ([2|) and are 

oo 

a s (uj) = q n A(uj - nu> m )a s [(uj - nw m ), 0] + 

n— — oo 

q n B{w - nuj m )a\ [(wj + nu m ), 0] 

oo 

a\(uji) = 2J r m C(uj — moj m )a s [(uj —mu} m ),0] + 



m— — oo 



r m D(uj - mu) m )a\\{wi + mu) m ), 0] 



With the modulators present, the spectrum at the 
signal and idler frequencies as would be observed as 
a function of position on the upper and lower screens 
of Fig (1) are S(u) = {2ir)(at(uj)a s (co)}, and = 
(2Tr)(al(uj)ai(ui)) , and evaluate to 
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q n \ 2 \B{uj - nuj m )Y 



r m \ 2 \C(uj + muj m )Y 



(4) 



The total number of counts at either frequency in 
gatewidth T is obtained by integrating Eq.((4|) over u. 

We take the coefficients of Eq.([T]) to correspond to a 
rectangular wavefunction with a temporal width equal to 
one unit 17f • These are A(u>) = D(cu) = 1, and B(u) = 
C*{ui) = exp(— ix) sm{x)/x, where x — (u> — luq)/2 and 
ujq is the center frequency of the signal spectrum. 

We assume ideal phase modulators with the functional 
form exp[ir5sin u) m t] at the signal and idler frequencies. 
The modulator coefficients in Eq. @ are then Bessel func- 
tions with q n = J n (-5 S ), and r m = J m (—Si). Fig. 2 
shows the spectrum of the signal frequency after passing 
through the modulator. Two limiting cases are of inter- 
est; in the first, Fig. 2 (a), the modulation frequency 
u> m = 0.1 is small as compared to the linewidth of the 
biphoton (1 unit), and the spectrum is nearly the same 
as without the modulator present. In the second case, 
Fig. 2 (b), the modulation frequency uj m = 10 is large 
as compared to the linewidth of the biphoton. In both 
cases we take the modulation depth 5 S — 2 radians. In all 
cases with ideal phase modulators the rate of generation 
of paired photons is R = l/{2n) J_ \B(to)\ 2 dco. 

To measure the two frequency correlation function one 
would observe a signal "click" on one screen, and rather 
immediately an idler click on the other screen. The po- 
sition of the idler click relative to the position u p — u s 
is recorded and integrated over all positions at the signal 
frequency. 

We take the two frequency correlation function, as 
observed in the frequency domain as (lj s , u>i) = 
(2tt) 2 (a\(u)i)al(uj s )a s (uj s )ai(uji)) . We expand by Wicks 
theorem and find that the result is the sum of two terms. 
The first term is the same as that obtained by a clas- 
sical (frequency domain) correlation of the intensities of 
the signal and idler. The second term is the result of 
time energy entanglement, and as described below ex- 
hibits what we term as nonlocal modulation. With a rel- 
ative detuning between the signal and the idler defined 



as A = [u 



(w. 



i)], and with T as the gate width, 



the classical term c(A) is 
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\B(u> — nuj m )\ \C[u) + A + mw OT ]| du> (5) 
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FIG. 3: Quantum correlation versus sideband number, (a) 
The signal and idler are modulated with the same phase so 
that 8s — Si = 2. Here the modulators act cumulatively, (b) 
The signal and idler are modulated with opposing phases so 
that 8s = 2 and Si — —2. Here, the modulators negate each 
other. In both portions of this Figure, w m = 0.1, as compared 
to the non-modulated biphoton bandwidth of unity. 

This classical correlation function is equal to the convolu- 
tion of the signal and idler spectra of Eq. ^ , is non-zero 
for all A, and is the same as that which would be ob- 
tained using any light source with the same spectrum. 

While the classical correlation function is continuous, 
the quantum portion is not; instead it is a set of delta 
functions that are located at integer multiples of the mod- 
ulation frequency u> m . With z = (A/u> m ), the quantum 
term q(A) evaluates to 



9(A) 
/(*) = 



T \ 

— J £ f(z)6[A + zu, m ] 
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FIG. 4: Here, the modulation frequency uj m = 10 is large as 
compared to the biphoton bandwidth of unity, (a) 8 S = 2 and 
Si = 2, and (b) S s = 2 and Si — —2. Cumulative modulation 
effects are no longer observed. 

(w TO = 10). Here, we find that the nonlocal modulation 
effects disappear. The requirement for cumulative modu- 
lation is that the modulation frequency be sufficiently low 
that all portions of the initial biphoton wavepacket expe- 
rience nearly the same phase. Since nonlocal modulation 
is the time-frequency counterpart of nonlocal dispersion 
compensation [21], the comparable requirement for non- 
local dispersion compensation is that the inverse of the 
rate of change of the dispersive constant with frequency 
be small as compared to the width of the frequency do- 
main correlation function. 

For ideal phase modulators, with R as the rate of gen- 
eration of paired photons [Eq.Q] we obtain the integrals 



g*r*_ n C(w -nu m )A*(u) -noj m ) 



du> 



If both modulators are absent, then Eq.(6) yields only 
a single delta function at zero. If a modulator is present 
in only one channel, the quantum term is a comb of delta 
functions. With modulators in both channels, and when 
the modulation frequency is small as compared to the 
spectral width of the biphoton wavepacket, the modula- 
tors interact in the sense that the cumulative modulation 
depth of the biphoton depends on the sum of the modu- 
lation depths of the independent modulators. 

Figure [3] shows the nonlocal correlation obtained with 
both modulators on. In Fig. [3] (a), the signal and idler 
channels are modulated with the same phase so that 
r5 s = Si = 2. The modulators act cumulatively, so as 
to produce the same correlation as would a single modu- 
lator with a modulation depth of four. In Fig. [3] (b), the 
modulators are driven with opposing phases, i.e. 5 S = 2 
and Si = — 2 . Here the modulators negate each other so 
that the effective biphoton modulation depth is zero. 

Of importance, the frequency of modulation in Fig. [3] 
is small (u) m = 0.1) compared to the bandwidth of the 
biphoton which is unity. Fig. 4 shows the correlation 
with the same conditions as Fig. [31 but this time with 



/oo 
c(A)dA = R 2 T 2 
-OO 

/OO 
q(A)dA = RT 
-OO 



(7) 



The first of these relations follows immediately by inte- 
gration of Eq.f5j). The second requires the assumption 
that |-A(o;)| is approximately unity, i.e., small parametric 
gain in the generating crystal.. From Eq.([7]), it is clear 
that to observe the quantum effects described here, that 
the gate width T be sufficiently small that RT << 1. 
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FIG. 5: Measurement of biphoton wave packets. Ampli- 
tude modulators in the signal and idler paths are driven syn- 
chronously. The coincidence count rate is measured by slow 
detectors and is plotted as a function of modulation frequency. 
The Fourier transform of this plot is the square of the mag- 
nitude of the biphoton wavefunction. 

We next consider synchronously driven amplitude 
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modulators and describe a Fourier technique for mea- 
suring biphoton wavepackets that are shorter than the 
temporal resolution of existing single photon counting 
modules. First note that because the incident biphotons 
arrive at the modulators at random times, the coinci- 
dence count rate will be unchanged unless the modula- 
tors are driven at the same frequency. We recognize that 
we may accomplish a time to frequency Fourier trans- 
form by multiplying in the time domain, i.e. by modu- 
lating, and integrating over time. In Fig. 5 signal and 
idler photons are amplitude modulated and are incident 
on detectors whose response time is long as compared 
to the temporal length of the photon wavepacket, and 
short as compared to the inverse pair generation rate. 
The modulation frequency at the signal and the idler is 
varied and the coincidence count rate is measured as a 
function of this frequency. With 5 s Si << 1, the inverse 
Fourier transform of the measured curve is the the square 
of the magnitude, or equivalently the quantum portion of 
the Glauber correlation function of the biphoton wave- 
function. 
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FIG. 6: Biphoton Wavefunction Measurement (a) Coinci- 
dence Counts per bandwidth versus modulation frequency 
with the average count rate set to zero, (b) Fourier Transform 
of part (a). 

The derivation is reasonably straightforward. We be- 
gin by multiplying the signal and idler operators at the 
output of the generating crystal by 1 + S s cos(w m i) and 
1 + SiCOs(u! m t), respectively. We form the correlation 
function = {a\(ti)a\(t2)a s (t2)ai{ti)) , Fourier trans- 
form and expand by Wicks theorem. Two non-zero terms 
are obtained: The first is the quantum term and the sec- 
ond is a classical term whose magnitude varies as the 
square of the generation rate. We take r = £2 — ti, 
average t\ over the modulation period, and integrate 
over r. With the biphoton wavefunction denoted by 
4>{lo) = A(uj)C* (w), and k = 5 s 8i/2-K, the quantum term 
to lowest order in k is 

/oo 
{<j)(uJ + uj m )(f)* (ui) + cc]dui (8) 
-00 

Eq. (J8]) is the inverse Fourier transform of the square 



of the magnitude of the biphoton wavefunction. Figure 
|6ja) shows the coincidence count rate as a function of 
the modulation frequency for a Gaussian biphoton wave- 
function with a duration of one unit and a delay of eight 
units in the signal channel. Fig. HJJb) shows the wave- 
function obtained by taking the inverse Fourier cosine 
transform of Fig. [Ha). We note that a different method 
for measuring short biphotons by correlation has been 
demonstrated by Silberberg and colleagues 
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